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Near a critical point, the equilibrium relaxation time of a system diverges and any change of 
control/thermodynamic parameters leads to non-equilibrium behavior. The Kibble-Zurek problem 
is to determine the dynamical evolution of the system parametrically close to its critical point 
when the change is parametrically slow. The non-equilibrium behavior in this limit is controlled 
entirely by the critical point and the details of the trajectory of the system in parameter space 
(the protocol) close to the critical point. Together, they define a universality class consisting of 
critical exponents — discussed in the seminal work by Kibble and Zurek — and scaling functions for 
physical quantities, which have not been discussed hitherto. In this article, we give an extended 
and pedagogical discussion of the universal content in the Kibble-Zurek problem. We formally 
define a scaling limit for physical quantities near classical and quantum transitions for different sets 
of protocols. We report computations of a few scaling functions in model Gaussian and large-A^ 
problems and prove their universality with respect to protocol choice. We also introduce a new 
protocol in which the critical point is approached asymptotically at late times with the system 
marginally out of equilibrium, wherein logarithmic violations to scaling and anomalous dimensions 
occur even in the simple Gaussian problem. 
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I. INTRODUCTION 

The study of critical points and their associated con- 
tinuum limits or field theories has been a central and 
enormously productive exercise in statistical mechanics 
and condensed matter physics. This study began with 
the problem of equilibrium finite temperature transitions, 
was then extended to their dynamics and thereafter to 
the interplay between criticality and finite size effects [T] . 
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This development has since been replayed in the theory of 
zero temperature quantum phase transitions [21 E]. Need- 
less to say, the development of powerful field theoretic 
methods, most notably conformal field theory in (1+1) 
dimensions [4^ ^.nd most recently the gauge gravity dual- 
ity [S]) have significantly advanced this line of work. In 
all these cases the field theoretic approach not only al- 
lows an understanding of the critical point but also the 
regions of the phase diagram adjacent to it. 

A new dimension to the study of the passage through 
critical points was introduced by Kibble |6j in the con- 
text of the expanding universe, since recast in the lan- 
guage of critical phenomena by Zurek [7| . Their proposal, 
the "Kibble-Zurek mechanism," is a theory of the defects 
generated in a system being cooled through a continuous 
symmetry-breaking phase transition at a small, but finite 
rate. The system inevitably goes out of equilibrium on 
the approach to the transition and arrives in the broken 
symmetry phase with different spatial regions realizing 
different orientations of the broken symmetry, and topo- 
logical defects as a result. The mechanism predicts the 
scaling of the number of these defects with quench time 
and has been successfully tested in liquid crystals [5]- 
[TU] , superfluids [TTtiTS] and superconductors (TU [T5] . It 
has also been recently generalized to the setting of quan- 
tum phase transitions by Dziarmaga |16| . Polkovnikov 
[T7j and Zurek et al. [TS] , where the role of temperature 
is played by a non-symmetry breaking control parame- 
ter, and to ramps across multi-critical points 19J and 
a gapless critical line [20] • The scaling of other physical 
quantities like excess heat with quench time has also been 
investigated [3T1. Polkovnikov and coworkers have also 
paid attention to the interplay between finite system size 
and parameter velocity in the context of quantum phase 
transitions. For reviews of this line of work, see [221 123| • 

The general problem posed by the work of Kibble and 
Zurek is that of a slow passage through a critical or multi- 
critical point which we shall term the Kibble-Zurek (KZ) 
problem (often referred to by the oxymoronic term "slow 
quench"). This problem is then characterized by a crit- 
ical point with its equilibrium physics and a "protocol," 
which is a particular path in the parameter space of the 
problem that touches the critical point. In the limit of 
asymptotically slow motion in parameter space, we ex- 
pect that the physics is dominated by the critical point 
and hence is, in an appropriate sense, universal. An im- 
portant task of theory is then to isolate this universal 
content and compute it. We should note that the study 
of universality in the KZ problem, especially in the quan- 
tum setting, is a part of a wider current study of non- 
equilibrium quantum dynamics; for broader perspective 
see the recent discussion of Polkovnikov et al. [21| . 

In previous studies of critical phenomena, two ideas 
have proven extremely useful. The first is the idea of 
the scaling limit in which various quantities of interest, 
such as thermodynamic densities and correlation func- 
tions, are postulated to obey certain homogeneity rela- 
tions. This set of scaling functions along with the critical 



exponents captures the full universal content associated 
with a given critical point. The second is the renormaliza- 
tion group which provides an understanding of the origin 
of this universality and a full computation of its content. 

In this paper we will make progress on the first front: 
We will formulate a scaling limit for the KZ problem and 
report model computations of the resulting scaling func- 
tions for a few classical and quantum problems. In appro- 
priate limits, these scaling functions will reduce to those 
in the equilibrium and the coarsening problems. Except 
for recent work by Deng et al. [25] and Biroli et al. [2^ 
discussing scaling functions in specific cases, the need 
for defining full scaling functions has been largely over- 
looked in the literature. Our contribution is to formalize 
the idea as a scaling limit for all physical quantities for 
any pairing of a critical point and a protocol. 

The second part of the program, beyond our ambition 
at present, would be the construction of a renormaliza- 
tion group flow. Such a program has been fruitfully pur- 
sued in sudden quench studies in classical models with 
stochastic dynamics (see 27J and references therein). We 
offer two modest steps in that direction here. First, in 
the classical context, we formulate the path- integral pre- 
viously written down only for sudden quenches to the KZ 
problem. This is, in principle, amenable to analysis by 
standard equilibrium renormalization group techniques. 
Second, we will prove universality with respect to pro- 
tocol choice for some model classical and quantum prob- 
lems. Specifically, we will show that the expectation that 
only the behavior of the protocols in the vicinity of the 
critical points is important is, in fact, correct. 

In the course of this paper we will offer a tripar- 
tite classification of possible protocols based on their 
"topology" — i.e. whether they cross, turn around at or 
end at the critical point. The scaling of quantities such 
as the defect density and the excess heat with quench 
time has been previously generalized [211 121] to arbitrary 
positive power-law behavior of the protocol on time near 
the critical point, that is, the first two kinds of protocols 
in our classification. Our contribution is the definition of 
a third class of power-law protocols that asymptotically 
end at the critical point and the identification of an in- 
teresting member — the marginal end critical protocol — 
that generates a one parameter family of non-equilibrium 
deformations of the equilibrium critical state worthy of 
further study. 

We turn now to the organization of this paper. We 
begin, in Section [ill by defining the protocols of interest 
and introducing the well known KZ time and length. In 
Section [Illj we lay out the scaling formalism for the KZ 
problem — the definition of the KZ scaling limit and the 
resulting scaling functions. Here we also discuss how the 
KZ scaling functions universally interpolate between the 
early equilibrium physics and the late time thermaliza- 
tion/coarsening physics. In Section IV we use the 0{N) 
vector model endowed with Model A dynamics as our 
classical model system and compute scaling functions for 
it in the Gaussian and large- iV approximations. We also 
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give an explicit proof of protocol universality for this 
problem. In Section |v] we turn to the quantum 0{N) 
model in the Gaussian approximation and compute var- 
ious quantities of interest here. We conclude with some 
remarks in Section VI. 

We note that a recent paper by Kolodrubetz et al. 
reports analogous results for the transverse field Ising 
model in (1+1) dimensions .29j . This work complements 
our own results on the quantum Ising universality class 
above the upper critical dimension (3 + 1) where the crit- 
ical theory is Gaussian iSlj. We comment on the one case 
untreated in |29j later in the paper. The present authors 
will also report computations of scaling forms in an in- 
teracting quantum field theory using the tools provided 
by the AdS/CFT correspondence elsewhere [30) . 



II. CLASSIFICATION OF PROTOCOLS AND 
THE KZ TIME AND LENGTH SCALES 

Consider a multicritical point in d spatial dimensions. 
Let {Oi} denote the set of relevant operators that cou- 
ple to conjugate fields {hi} and have scaling dimensions 
{Ai}. At the critical point, we set (Oi)eq = 0. The set 
of scaling dimensions {A.^} and a host of scaling func- 
tions constitute the static universal equilibrium content 
of this critical point. Along with the dynamical exponent 
z obtained from the time-dynamics, they determine the 
universality class of the critical point. Let S parameterize 
a path in the space of conjugate fields {hi} such that S is 
zero at the critical point. For every such path, define the 
correlation length exponent ly to be l/(d— A), where A is 
the scaling dimension of the most relevant operator that 
has a projection along the path near (5 = 0. Along the 
path, the correlation length then diverges as 6^'^ close to 
the critical point. 

A useful example to keep in mind, particularly for the 
next section, is the ferromagnetic critical point. The rel- 
evant operators are the scalar energy operator and the 
vector magnetization, coupling to the conjugate fields of 
temperature and magnetic field respectively. 

Unless mentioned otherwise, all length/time-scales are 
dimensionless and are measured in units of some micro- 
scopic length/time. We also set h and fc^ to 1. 



A. Protocols in 5 

Consider a system prepared in equilibrium at t = — c» 
at some fixed distance 5_ > away from the critical 
point being evolved in time along the path S{t) in con- 
jugate field space. For simplicity, we restrict ourselves 
to paths with a unique tangent most relevant operator 
near the critical point. The symmetry and the leading 
order behavior of S{t) in time near the critical point then 
completely classify the different universality classes that 
arise in KZ dynamics. The symmetry of the path is de- 
termined by the symmetries of the most relevant opera- 



tor: if it respects all the symmetries of the critical the- 
ory, the path is non-symmetry breaking (NSB), else it 
is symmetry-breaking (SB). Below, we carefully classify 
the leading order behavior of 6{t) into three categories 
depending on whether the power of t near 5 = is odd, 
even or negative. 

• Trans-Critical Protocols (TCPs): These protocols 
take the system across the critical point. They 
smoothly interpolate S between S- > as i — ^ — oo 
and (5-)_ < as i — > cxi, crossing the critical point at 
i = 0. An example of a TCP is 



-Sq tanh ■ 



t 



(1) 



We will shortly show that the dynamic scaling func- 
tions are universal with respect to the behavior of 
the protocol near t = 0. In anticipation of this re- 
sult, we classify the entire family of analytic TCPs 
by their leading order behavior in a time-scale r 
near i = as 



Sit;T) 



6- 




t — > — OO 
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t — > oo 



(2) 



with a odd. a = 1 is the linear protocol, = 3 the 
cubic and so on. 

• Cis-Critical Protocols (CCPs): These protocols 
keep the system in a single phase and touch the 
critical point aX t — Q. They smoothly interpolate 
5 between 5^ as t ^ — oo and > as i — >■ oo 
through the critical point at t = 0. An example of 
a CCP is 



5{t; t) = 5q tanh^ 



t 



(3) 



Eq. ^ with a even classifies CCPs. a = 2 is the 
quadratic protocol, a = 4 the quartic and so on. 

• End-Critical Protocols (ECPs): These protocols 
keep the system in a single phase and asymptot- 
ically approach the critical point. They smoothly 
interpolate 5 between (5_ at i — )> — oo and zero as 
t oo. An example of an ECP is 



m 



(5_(e 



-t/T 



1) 



t/r + e-*/^ 



(4) 



ECPs may asymptotically approach the critical 
point with or without a time scale (exponential vs 
algebraic). We restrict ourselves to the family of 
scale-free protocols: 



m 



\a\ 



t — > — OO 

t > r , 



(5) 



where a < and r is a microscopic time-scale sep- 
arating the early and late-time behavior of 5{t). 
Unlike in the above two cases, a is not required to 
be integer valued. 
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FIG. 1: From left to right, examples of Trans-Critical (TCP), Cis-Critical (CCP) and End-Critical (ECP) protocols. The 
critical point is at t = 0. The leading order expansion of 5 is valid in the shaded region. 



B. The KZ length and time 

A system evolving from t — — oo by a TCP or a CCP 
must fall out of equilibrium near the CP due to critical 
slowing down. This is signaled by a diverging time scale, 
^t, in unequal time equilibrium correlation functions, 
is interpreted as the relaxation time near a classical CP 
and the inverse local gap to excitations near a quantum 
CP. The time at which the system falls out of equilib- 
rium is defined to be the KZ time t^- The KZ time 
defines a KZ length Ij^ t]/^ . Ij^ manifests in dynamic 
correlation functions as a crossover scale between equilib- 
rium and non-equilibrium correlations. More intuitively 
[6l[22], in an ordering transition, order is unable to form 
on scales larger than due to the finite quench rate, and 
domains of broken-symmetry phase of size Ik persist in 
the ordered phase. We now separately consider the cases 
of TCP/CCP and ECP. 

TCPs/CCPs. The scaling of the KZ length and time 
with T follows essentially from considerations of [T51 155] , 
which we recapitulate here. A purely quantum version 
can also be formulated using Landau-Zener arguments 
and perturbation theory [T71 [TH]. We define ^{t;T) and 
(t; t) as the equilibrium correlation length and time at 
S{t;T). When t <^ — r, S{t;T) decreases slowly from 
and the system evolution is adiabatic/quasi-static. De- 
fine ts to be the time-scale over which <5(t;r) changes 
appreciably in the window r near t = 0. Using Eq. 
we find that = (5/|(5| ^ ts therefore decreases as 
t — > 0. The time scale ^t, on the other hand, diverges as 
\5o{t/T)'^\~'^^ . At t = —Ik, the two scales cross and the 
system must fall out of equilibrium. This last considera- 
tion fixes Ik'- 

UtK]T)^-tK, (6) 

which implies 




ECPs. The non-equilibrium behavior here is not 
signaled by the KZ time. Instead, the scale-invariant 



ECPs defined in Eq. ([5| fall into three classes depend- 
ing on whether \a\ is greater than, equal to, or lesser 
than \ac\ = \livz). When \a\ > \ac\, > ts for all 
late times, with the inequality getting stronger as S ap- 
proaches its critical value. The system is far out of equi- 
librium and its physics is not controlled by the critical 
point as t — > oo. Similarly, when \a\ < jad, the sys- 
tem is arbitrarily close to equilibrium sufficiently near 
the critical point. \a\ = jad is the marginal case when 
the instantaneous relaxation time of the system and ts 
differ by a factor which becomes constant at late times. 
£,tit) ~ t and £,{t) ~ t^/^ should then play the role of a 
KZ time and length in the scaling forms. 



III. THE KZ SCALING LIMIT 

The introduction of t^ and Ik is reminiscent of the 
introduction of finite size cutoffs in the theory of equilib- 
rium critical behavior, and we are led to analogs of the 
finite size scaling limit and finite size scaling functions. 
We define the KZ scaling limit to be the limit t — > cxo 
when time and length scales are measured in units of 
the diverging KZ scales, and Ik- S{t;T) — >• in this 
limit, and the system is arbitrarily close to the CP, evolv- 
ing non-adiabatically for all t/i^ (the shaded region in 
Fig. ^p) . In the case of the marginal ECP, the relevant 
length and time scales at some late time t are the instan- 
taneous correlation length and time scales, i.e. ^(t) and 
£_t{t) respectively. The scaling limit is then defined as the 
limit t oo holding x/^(i) fixed. 

We now turn to the definition of the scaling functions 
that arise in the KZ scaling limit. We discuss the scaling 
functions of various physical quantities and their asymp- 
totic forms in specific cases. Any readers new to scaling 
theory are encouraged to arm themselves with the Model 
A example in the Gaussian limit discussed in Sec. |IV| 
and interpret the definitions below in that context [T]. 
We initially consider classical critical points and discuss 
asymptotic forms when the transition is ordering. We 
then comment on the special features of quantum criti- 
cal points. 

Our notation is to denote the absolute value of any 
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vector k by fc, and to use x and i respectively to refer to 
the scaled length and time, x/Ik and t/tj^. We reserve 
calligraphic lettering for scaling functions. We also as- 
sume that in equilibrium, the symmetry is spontaneously 
broken across the transition in the NSB direction. 
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FIG. 2: Left: 5 as a function of t close to t = with KZ 
times shown for a slow (2, dashed) and fast (1, solid) quench. 
In the slower quench, the system falls out of equilibrium ear- 
lier {\tK2\ > I^Kil) but at a smaller distance from the CP 
{5{—tK2) < 5(— ^Ki)). Right: The correlation length of a 
system quenched at a finite rate (soHd) vs ^(i; t) (dashed). 
KZ split dynamics into (i) Adiabatic (ii) Sudden and (iii) 
Post-quench (here coarsening) regimes. The KZ scaling limit 
describes (ii) with (i) and (iii) as asymptotes. 



straightforward extensions of the form in Eq. ([8]). Note 
that Qo can be identically zero. This is the case, for in- 
stance, in a zero-field temperature quench through the 
ferromagnetic critical point when O is the spin operator. 
The finite quench rate prevents order from forming on 
scales longer than the KZ length, and (0)t, the average 
magnetization, remains zero at all times. 



2. Asymptotic form near equilibrium 

By construction of the protocol, we should recover 
equilibrium scaling forms in certain limits as t ^ ±oo. 
However, the precise limits are subtle and we derive them 
below. Recall that the equilibrium scaling limit is the 
limit of (5 — holding x/^ and {t — t')/^^ fixed, wherein 
the two-point correlator has the scaling form, 



(o(x,t)o(o,<'))r 



? yoo 



X t-t' 



(10) 



The ~ symbol here is distinct from that in the previous 
subsection. 

We consider the limit in which the KZ scaling form 
Eq. ^ reduces to the equilibrium scahng form Eq. (10 1. 
First, we observe that the relation. 



x/Ik C(i;T) 



x/at;T) 



t 



(11) 



A. TCPs and CCPs 

1. Scaling forms of correlation functions 

Consider the scalar operator O with scaling dimension 
A. We assume that the theory has translational and 
rotational invariance [52]. The KZ scaling forms for the 
one and two-point connected correlation functions are 



{0{^,t))r = Go{t;T) 



-yo 



/A - - \f 

l-K \ ''K 

(0(x, t)0(x', t'))r = Goo(|x - x'l, t, t'; t) 

^ p^yoo 



I,. 



(8) 



When we use ^ to indicate a scaling form of a correlation, 
as in (|8|, we have a precise limiting statement in mind. 
For example, in the case of the two-point connected cor- 
relator, what we mean is 



lim fJ^Goo (x, t, t';T) G. 

>-*oo 



X t t' 
oo \ -j—, —, — 

f'K 



(9) 



Ik ' 



where r — )■* oo means the limit where r — > (X) with 
and 5^ held fixed. 

lit 

The scaling forms of all higher order cumulants 
and cross-correlators with other relevant operators are 



implies that x/^ is fixed whenever x and t are held fixed. 
Similarly, t/^^ and i'/C^ ^-Iso fixed by t and t' . Thus, 
an alternative to the KZ scaling form Eq. ^ involving 
the same two arguments as the equilibrium scaling form 
may be given: 



lim ^ 

"— f * oo 



,{x,t,t';T) 



C(2) 
Voo 



X t-t' t + t' 



t. 



(12) 



Physically, when the correlation time S,t is smaller than 
the time-scale of change in S, i.e. \t\ ^ t^, the system is 
in instantaneous equilibrium on length and time-scales, ^ 

(2) 

and respectively. This is the limit in which Qqq must 
reduce to Goo- Thus, in the limit {t + 1') ±oo holding 
x/C and {t-t')/^^ fixed: 



yoo 



X t-t' t + t' 



yoo 



X t-t' 



(13) 



For the original KZ scaling form Eq. (|9|, this translates 
to the requirement. 

Goo i, i') ^ P""-^ Gz (x (i - i') i^'^n . (14) 

in the same limit. As expected, £.{t;T)/^{t';T) — >• 1 so 
that there is a single diverging length in the equilibrium 
system. Furthermore, time-translation invariance is re- 
covered. 
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In the example of a temperature quench through a 
ferromagnetic critical point, the equilibrium equal-time 
connected spin-spin correlation function decays exponen- 
tially on a length scale ^ on either side of the critical 
point; consequently, Eq. (14) when t = t' must asymp- 
tote to 



P""^ exp( 



3. Asymptotic form for coarsening dynamics 

It is generally believed that a system quenched to an 
ordered phase with multiple vacua undergoes coarsening 
|53j , whereby each local broken-symmetry region grows in 
time and the system is asymptotically statistically self- 
similar on a characteristic length scale, lco{t)- Put an- 
other way, the two-point function heals to its equilibrium 
value on the scale ^ within each "domain" , and is expo- 
nentially suppressed between domains, each of growing 
length ^co S> ^. In the late time regime, dynamical scal- 
ing is expected to hold when there are no growing scales 
competing with /co- For example, 



(o(x,t)o(o,t')>5^e 



-2Ar'co 



. ''CO 

as t,t' — )■ oo. Sudden quench studies in the Glauber- 
Ising model the 0(A^)-model with N = oo [32] and 
numerical simulations in other systems have verified this 
hypothesis. 

When present in the KZ problem, we expect coarsen- 
ing physics to emerge deep in the ordered phase, i.e. as 
t,t' — > oo, on the length scales lco{t) , lco{t') ■ Proceed- 
ing as in the previous sub-section, we conclude that as 
t, i' —i' oo holding x/l^oit) and x/lco{t') fixed, the two- 
point function must have the limiting form: 



Gooix,i,i') 



where 9 



~2aiyA /^co 



avz + 1 

Zd 



(15) 



Zd is a dynamic exponent specific to coarsening and we 



have used the relation, 



U(t;r) _ 
Ik ~ 



Note that the 



limiting equilibrium form requires holding x/(,(t) fixed 
and that lco{t)l£,{t) diverges as t — > oo. 

In the simplest case of Model A dynamics for spins 
with A'^-components, Zd = z = 2 and lco{t) ~ Vi in the 
infinite- limit. Although this phenomenology holds for 
systems with and without topological defects, the specific 
form of the equal-time spin-spin correlation function de- 
pends on the presence of topological defects; for more 
details, see the excellent review [551 . 



4- Scaling form of the non- equilibrium correlation length 

The non-equilibrium correlation length is defined 
to be the inverse of the decay constant on the longest 



length scales of the two-point equal-time correlator in real 
space. In equilibrium, this length is the cross-over scale 
in correlation functions between fluctuations dominated 
by one fixed point and another. A particularly simple 
definition of ^ne is. 



Cnc(i;T) = 



/ / dx<^x'^Goo{x,t,t]T) 
J dx'^Goo{x.,t,t]T) 



(16) 



and is useful when Goo is always positive. A more gen- 
eral definition is through the smallest imaginary part of 
the poles of Goo{k,t,t;T) in fc-space; we will elaborate 
on this when the necessity arises. Biroli and coauthors 
[25] recently discussed the KZ scaling form of ^ne in NSB 
TCPs. They observed that ^nc must asymptote to the 
equilibrium correlation length, S,{t]T), as i — ^ — oo; that 
it must scale according to the critical coarsening form, 

l/z 

til J when |i| ~ 0(1); and that it must asymptote to the 
coarsening length ^co as i — > oo. Their proposed scal- 
ing form can be derived from that of Goo and can be 
re-written as 



Cnc(^; '^) ^ ^K^n 



The asymptotic behavior of ^nc for i 
translates to the limits 



(17) 



±(X) formally 



\x\' 



X 
X 



—oo 
oo. 



(18) 



Absent coarsening physics, the right asymptote repro- 
duces the instantaneous correlation length. 



5. Scaling form of the number of defects 

An intuitively appealing picture of the lack of order on 
length scales greater than the KZ length is through topo- 
logically protected point defects of characteristic separa- 
tion Z^, and/or defects of dimension p with characteris- 
tic separation 1^ in any hyperplane of co-dimension p. 
This picture is really only meaningful if the separation 
between defects is much larger than the equilibrium cor- 
relation length; a constraint met only in the coarsening 
regime, t ^ t^. In the coarsening regime, the density of 
a defect of dimension p in the hyperplane of co-dimension 
p should scale as l/lco{t;Ty~^ ■ At fixed positive i, this 
reproduces the celebrated scaling of the defect density 
with r, 



Density of defects 



d—p 



r^^^. (19) 



The above scaling with r has been verified in experiments 
in liquid crystal [10], superconducting [T5] and ion-trap 
[31] systems. Although we have presented the above 
as a natural scaling ansatz, it can be derived from the 
known scaling of correlation functions, for example, by 
the methods of Halperin-Liu-Mazenko for classical tran- 
sitions [351 IMI- 
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6. Scaling forms of thermodynamic analogs 

The dynamics of classical systems is typically modeled 
phenomenologically by stochastic differential equations, 
possibly with conservation laws [37] . ft is often useful to 
reformulate the stochastic dynamics in d-dimensions in 
terms of a path integral in (c?+ 1) [3Hlini]- For Gaussian 
noise, the generating functional of correlation functions 
of the fundamental field 0{x,t) is 



Z[J] 



'DO{x, t) exp 



(fxdt{C + J{x,t)0{x,t)) 



Goo(.x,t,t') 



X t 



(24) 



where now the precise meaning of the asymptotic equality 
is 



lim 

, A fixed 



t-Goo{x,t,t')^gc 



X t 

t^' t' 



We emphasize that we cannot appeal to results from equi- 
librium theory or sudden quenches to/near the CP to 
understand the asymptotic forms of correlators here. 



Here J is the source for O, and £ is a local Lagrangian 
density in which parameters of the protocol (like r) ap- 
pear as couplings. When hyper-scaling is obeyed, all KZ 
scaling forms follow from the scaling of the associated 
free energy. 

In this formulation, it is natural to define a time- 
dependent free energy density. We divide the (d + 1)- 
dimensional space-time into a stack of spatial slices with 
volume L'^ and temporal length At. Neglecting boundary 
effects, we compute the free-energy density of a slice: 



/(t;r)= hm 



-log(^) 



(20) 



where 



Z{t) = [ VOexp i-f d'^x ( 



t+At 



dtC 



For time-independent 6, f is also time-independent, and 
its leading non-analytic dependence on S is /na ~ 6"'^. 
When S varies with time, we conjecture a scaling form 
for fna.(t', t) in the KZ limit : 



/na(i;r)^-^(- 



(21) 



When t/tjc <C — 1, the system is asymptotically in equi- 
librium. Thus, 



when t/tK (22) 



This is also the expected asymptotic form of J- when 
t/iK > 1 for a SB quench or a NSB CCP. 



B. Marginal ECPs 

The roles of the KZ time and length are played by 
it{t) ^ t and i{t) - for the marginal ECP. The 
scaling forms of the 1-point and 2-point unequal time 
correlation function are 



Gait) 



1 



(23) 



C. Quantum systems 

We comment below on the major differences that arise 
in the treatment of quantum systems. 

• The scaling of the correlation functions proceeds 
as before. However, the analogs of thermodynamic 
quantities of interest for protocols that begin with 
the system at T = are now the excitation en- 
ergy density in excess of the energy density in the 
adiabatic ground state, or "heat" density [21] [54] . 
and the entropy density for which plausible defini- 
tions can be constructed from the diagonal entropy 
[40l |4T] or the entanglement entropy of a macro- 
scopic subregion. These are respectively expected 
to exhibit the scaling forms 



q{t\T) 

S{t-T) 



1 



id+z 
1 



t 



l<i~^ \ t 



(25) 
(26) 



For an isolated quantum system, s and q are con- 
stant whenever the evolution is adiabatic, i.e. as 
t — )■ ±oo provided the vacuum is unique. When 
the vacuum is not unique as i — cx) and the system 
exhibits coarsening, the evolution is weakly nonadi- 
abatic in that most of the degrees of freedom evolve 
adiabatically. Evolution by the marginal ECP is 
not adiabatic by construction. 

• For protocols that begin with the system in equi- 
librium at T > 0, a new dimensionless parame- 
ter, ksT/^h/ti^), now enters the quantum problem. 
Along with the quantities held fixed as the quench 
time is taken to be arbitrarily large in Sec. |III A 1[ 
we also hold Tt^ fixed. For example, the scaling 
form for the 2-point unequal time correlation func- 
tion is now 



(O(0,<)O(x',i'))r,T 



1 



X 



The definition of the entropy goes through as before 
but the excitation energy density is now measured 
with respect to the energy density that would be 
obtained in a strictly adiabatic evolution [55) . 
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Integrable quantum systems allow the existence of 
sharply defined quasiparticle excitations; their den- 
sity resolved by momentum can serve as a uniquely 
quantum observable for CCPs and ECPs. For NSB 
TCPs, it is not straightforward to define quasipar- 
ticles and separate them from domain walls and 
other topological defects. Currently, we are not 
aware of an integrable system above one dimension 
where this question can be properly posed (in 
one dimension there is no real distinction between 
quasiparticles and domain walls). Note that the 
definition of the thermodynamic quantities of heat 
and entropy do not rely on integrability. For in- 
stance, in a cyclic process like a CCP, q{t; r) when 
< ^ t/r is the difference of the system's energy den- 
sity between symmetric time points. 



D. Comments 

We conclude this section with two comments. 

• The Kibble-Zurek picture was initially proposed for 
the NSB-TCP through a thermal transition. The 
time-dynamics was decoupled into three regions: 
the adiabatic regime for t < —t,^, the diabatic or 
sudden or impulse region from —Ik to when the 
system is frozen, and the post-quench regime for 
t > when the system is unfrozen and evolves 
through domain growth [55] , defect-anti-defect an- 
nihilation etc. Recently this picture was ex- 
tended to account for evolution in the impulse 
regime using critical coarsening results. The KZ 
scaling forms introduced in this article also probe 
dynamics in the impulse regime with the adiabatic 
and the post-quench regimes acting as asymptotic 
limits. 

• A finite system dimension L can be readily accom- 
modated in our scaling forms in the combination 
L/Ik- The results obtained through adiabatic per- 
turbation theory [ITl |43] can therefore be fitted 
within this framework. 



IV. CLASSICAL SYSTEMS WITH MODEL A 
DYNAMICS 



dt 



OF 



(27) 



X and t are dimensionless and measured in units of the in- 
verse cutoff and A~^ respectively. is a zero-mean 
spatially uncorrelated white-noise stochastic variable for 
every a. The variance of (a is chosen to be 2 so that the 
long-time limit of the structure factor computed from the 
equation of motion when F is time-independent is equal 
to the equilibrium structure factor: 

(C„(x, t)CMx', t')) = 26^pS^{^ - ^')d{t - t') . (28) 



A. Gaussian Limit 

In this limit, we drop the (j)'^ term in F. The critical 
point is at rg = 0, /i" = 0, i.e. at the origin in the (ro, h") 
parameter space. Let us restrict ourselves to paths in this 
space that lie along the axes and include the origin. The 
equilibrium theory is sensible only at the origin and when 
To > 0. We can therefore study CCPs and ECPs along 
the non- negative vq axis {S = rg) ■ Luckily, we can also 
study TCPs along the same axis as the time-dependent 
fields are finite even when rg < 0. This physically un- 
interesting protocol is pedagogically useful. The critical 
exponents for all such paths are f = 1/2, z — 2. In the 
remainder of the discussion, S = ro and h°',u — 0. 

Eq. (27) is linear in (j) and hence diagonal in Fourier 



space. As it is a first-order equation, it can be solved for 
any protocol S{t;T) = ro(i;r), a given initial configura- 
tion at i = — oo and a fixed noise realization: 

<^(k,t) -0(k,-oo)+ r dt'e-/*'^*"('='+'-«(*"'"»C(k,i') 

J — OO 

(29) 

We have dropped the component label of (j) for brevity. 
The equal time structure factor for each component, de- 
fined as (0(k,t)0(k',t)) = {2TT)^S^{k + k')G^^{k,t;T), is 

+ 2 f dt'e-2/,U*"(fc^+n,(t";r))_ (30) 
G^^{k, — oo) is assumed finite for all k. 



We illustrate the universality of the KZ scaling limit 
and explicitly compute scaling functions of a vector op- 
erator (f) in the simplest setting: Model A [37] dynam- 
ics with a Landau-Ginzberg- Wilson (LGW) free energy. 
Model A dynamics is dissipative and obeys no conserva- 
tion laws. Let be an iV-component vector field in d = 3 
dimensions. The dimensionless LGW free-energy and the 
equation of motion are, respectively. 



F = / d^x 



1 



1. Scaling limit for tanh TCP 

The universality of the scaling limit with respect to 
details of the the protocol is already apparent when we 
consider the simple TCP ro{t;T) = — tanh(f/T). This is 
an example of a protocol that is linear near the critical 
point, i.e. a = 1 in Eq. ([2]). Consequently, the KZ time 
and length scales are 



.1/4 
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Re-writing Eq. (30) in units k = kl^ and i = t/tf 



G^^{k,t;T) = G^^{k/lK, ~oo) 



+ 2t. 



As T — >■ oo, tanh{t" /t^) — > t" and the KZ scaling form 
of the two-point equal time correlation function is 



G^^{k,t;T) - llg^,^{k,i) 



ll / 2df'e-2/<^'^*"(^'-*") 



(31) 



Observe that the scaling function only depends on the 
leading behavior of the protocol near the CP. For all pro- 
tocols such that ro(t;T) ^ ^i/'T hi the vicinity of zero 
time, the scaling form of the structure factor is the ex- 
pression that we have just derived. 



2. Scaling functions for all TCPs and CCPs 































I 






1 







FIG. 3: G4,4,(k,t) vs k at fixed time slices for the linear 
TCP. The blue, red and green solid lines are respectively at 
t — 0.2,0 and —0.2. The green dashed line is the correlator 
if the system were in equilibrium at t = —0.2. Inset: The 
scaling form of the non-equilibrium correlation length vs t in 
solid. The dashed line is the instantaneous correlation length 
in units of Ik, £,{'t\T)/lK- 



Wick's theorem informs us that all higher order cumu- 
lants of 4> only depend on G^^ . The details of the proto- 
col enter the expression of G^^ only in the combination 
t^TQlt] t). In the scaling limit, for any a in Eq. ([2]), it is 
easily seen that 



(32) 



The scaling function of the equal time correlation func- 
tion for any protocol with the near-zero behavior in 
Eq. ([2| is therefore 



Q,,{k,i) = 2 



^ ^2fi,di"{p+{-i"r) _ 



As seen from the Fig. [3] for the linear TCP, G^^ and 
£no match equilibrium forms as < — )• — oo and A; 3> 1: 



(-t)i/2 



(33) 



Observe that Q^^ is finite at A; = when the system is 
near and exactly at the CP, indicating the suppression 
of order on length scales longer than l^- As £ — oo, the 
non-equilibrium correlation length grows without bound 
for TCPs due to the inverted potential. For CCPs, on 
the other hand, all scaling forms asymptote to the equi- 
librium forms at large positive times. 

In Appendix |A] we compute the partition function in 
(3-1-1) dimensions and demonstrate the validity of the 
scaling hypothesis for the free-energy density /na- 



3. Scaling functions for the marginal ECP 

The protocol ro(i) = rt^jt defines the marginal ECP 
at late times for the Gaussian theory. The scaling limit 
is independent of the short-time regularization of 7'o(t) 
and the initial field configuration. The asymptotic form 
of G^^ in the Gaussian limit is 

G^,(k,t)^tg,^(,k^ft) 



2a „4ro + l 



dye ^ y 



The marginality of the deformation is manifested in the 
distinction between ^/^^(x) and Q'^'^^x) when we make the 
identification ^(t) = \/t/rQ. Furthermore, at late times, 
the non-equilibrium correlation length approaches a 
finite multiple of the instantaneous correlation length: 



TQ + l 



B. Large- TV limit 

In the infinite N limit, the LGW theory is exactly solv- 
able in terms of 1 and 2-point correlators. All higher 
order correlators follow by the application of Wick's the- 
orem. The 1 and 2 point correlators are known exactly 
in equilibrium and can be reduced to quadrature with 
Model A dynamics. The theory in d = 3 that we discuss 
here exhibits non-trivial critical behavior and affords us 
a probe of SB paths and coarsening physics. 

Assume in Eq. (27) that /i* ^ for i = 1 and zero 
otherwise. This specifies all one-point correlators for 
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I > 1 to be zero. The two-point correlator for every 
component is the same and is a function of the variable 
= vq + u{(p'^). The well-known self-consistency equa- 
tions relating G^^ and rri^ are 




2 G^^ (fc) 



N 



(34) 



(35) 



where A is a cutoff on the maximum allowed |fc|. 

The third equation that completes the theory in equi- 
librium is, 



fc2 



(36) 



The critical point in the equilibrium theory is at tq — 
rc — —uA/27t'^ and — 0. The two relevant operators 
that couple to vq and ft." are respectively and 0ct, 
with scaling dimensions 1 and 5/2 in d = 3. They are 
respectively NSB and SB. Consequently, we can study all 
three kinds of protocols along the vq axis {i' = 1, z — 2) 
or along any of the h" axes {i^ = 2/5, z = 2). 

The equation of motion ( 27 ) for Model A dynamics is 



,(k,t) 



dt 



-{k' + m')Mk,t) + ft" (2^)^,5^ (k) + C 



(37) 



This specifies the functional dependence of G^^{k,t]T) 
on and completes the dynamical theory. 



Scaling limit for NSB TCPs and CCPs 



One can now re-express ( 34 ) as 



G^^{k, t; r; 0) 



d^k 



(27r)3 



G^^{k,t;T;m^) 



(40) 



Let us now take the scaling limit. The KZ length 

and time scales are = (t'^/Jq) and Ik — \ftK- 
Dimensional considerations imply that the scaling form 
of 'm?{t\T) is m'^{t;T) ~ M'^{i)/P^. Using the leading 
order behavior of the protocol near the critical point, 
S{t; t) ^ {—iy/lK^ we derive the following scaled form of 



Eq. (40) 



{~iY + 2u / di' 



(27r)3 



-2ji, di"{P+M^) 



We note the following points: 

1. The cutoff A appears in the unsealed equation only 
as the upper limit on the momentum integral. In 
the scaling limit, the upper limit ^ Al^ — >■ oo and 
the scaled relation is cutoff-independent. 

2. The scaled relation, being a limiting form, is sim- 
pler than the unsealed one. 

3. The relation above is unchanged by keeping higher 
order terms in r'o(t; r). This is tantamount to prov- 
ing the universality of all correlators with respect to 
details of the protocol away from the critical point 
in the scaling limit. 



For NSB paths, ft" = and 6{t;T) = ro(t;r) - r^. 
Eq. ( 35 ) implies that the one-point correlator for all com- 
ponents is zero: = 0. We henceforth drop compo- 
nent subscripts. We will make frequent use of the Gaus- 
sian result (30), which we reproduce here for the reader's 
convenience: 

G^^(fc, t; t; ro) = G^^{k, -oo) 

dt"(fc2+ro(t";r)) 



The notation G^^{k,t;T;ro) emphasizes that G^^ de- 
pends on ro(t; r). 



Observe that the solution to Eq. ( 37 ) is exactly that of 



the free Gaussian case — Eq. (29) — with ro(i;T) replaced 
by 171^ {t; t) . Consequently, the third equation that com- 
pletes the dynamical theory is 

G^^{k,t;T) = G^^{k,t;T;m^) , (38) 

where the right hand side refers to the Green's function 
of (38) with ro(t;T) replaced by m?{t\T). The critical 



value Tc of To in this notation is 

d^k 



{2n) 



G^^(k, t; t; 0) . 



(39) 



2. Scaling forms for the linear TCP 



Fortunately, Eq. (41) can be solved when a — 1. The 



choice It = v87r simplifies pre-factors and makes appar- 
ent the form of the solution. Defining 

fit) ^ e^P.^i'^^i') or M\i) = (42) 

the solution is 

/(£) = -3i/3r(l/3)e2*'/3(t Ai(P) + Ai'(P)) . 

Ai is the Airy function of the first kind and F is the 
gamma function. We may reduce Q^^ to quadrature: 

Cy,,(fc,f) = 2 r dt'e-^/^-i^-C^^+AI^ (*■")), 



and compute all higher cumulants using Wick's theorem. 

As t — >■ — cx), we recover the equilibrium behavior in 
Eq. (33) with ^ ~ l/(— t). The results from the coarsen- 
ing literature [33] in this theory are : 
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FIG. 4: G4,<t,{k,t) vs k at fixed time slices for the finear TCP. 
The blue, red and green solid lines are respectively at t = 
0.5,0 and —0.5. The green dashed line is the correlator if the 
system were in equilibrium ait — —0.5. Inset: jCno vs t (solid) 
and £,/Ik (dashed). 



We should therefore expect that Cne ~ and Q^^ ^ 
i^l"^ exp(— 2fc^t) as i — >■ oo. As is seen from Fig. [ij both 
asymptotes are correctly predicted for £nc- We have 
checked this for Q^^ as well. This verifies that the asymp- 
totes predicted in Eq. ( 18 ) and in Eq. ( 15 ) are correct in 
this theory. 




3. Scaling limit for SB TCPs and CCPs 

Setting ro = rc, we can explore symmetry-breaking 
protocols along the axis. The implicit relation from 
the self-consistency equations for the linear TCP is 



,iJJtzm, (43) 

|t-i'|3/2 



Although scaling is guaranteed, we cannot proceed fur- 
ther as the solution to this equation is not known. 



V. GAUSSIAN QUANTUM FIELD THEORIES 

We now turn to the KZ problem near quantum critical 
points. Again, it will prove instructive to investigate the 
class of ferromagnetic critical points with 0{N) symme- 
try. We will consider the simplest non-trivial case, that 
of Gaussian scalar fields. The Lagrangian we consider is. 



(44) 



where we permit to depend on time {6 = m^) . Stan- 
dard critical exponents include z = 1 (owing to relativis- 
tic symmetry when m is a constant) and v = 1/2 (owing 



to the Yukawa form of the Green's function for static 
sources when m is a constant). 



A. Second quantization 



The second quantized treatment of 
expansion 

d'^k 



is based on the 



(x,t) 



(27r)" 



where we set 

0k(^) — Jk\''J^k V J-kV-y^-k 

and impose the commutation relations 



(45) 



(46) 



(47) 



Owing to rotational symmetry, the mode functions fk{t) 
and Green's functions only depend on the magnitude of 
the momentum k. The mode functions fk{t) satisfy the 
mode equation 



fk{t) = 0, 



where 



nl{t) = fc2 + m2(i). 



One must also impose the Wronskian condition 

fkfk ^ fkf'k = * 



(48) 



(49) 



(50) 



in order to obtain standard commutation relations be- 
tween 0k and its conjugate TTk = (j)]^- Once the mode 
functions fk are specified, a Fock space vacuum |0) can 
be defined through the conditions ak|0) = 0. 

In general, Eq. (48 1 is hard to solve exactly. However, 



the techniques of WKB provide an approximate solution 
when flk is real, positive and varying "slowly enough"; to 
this end, let us define the oscillatory, positive frequency 
WKB solution: 



.hit) 



1 



exp 



i / dt'^kit') 



(51) 



where the lower limit of integration can be specified 
at our later convenience. If one defines dk(t) through 
the equations <j)k{t) — fk{t)ak{t) + f^{t)a'_^{t), and re- 
quires [di,{t),al^{t)] = (27r)''(5''(k - k') for aU t, then 

if kit),. f kit)) and ifkit)Jlit)) are related through the 
standard Bogoliubov transformation. 

As is slowly varying and positive at large negative 
times, the evolution is adiabatic and we require that the 
fk coincide asymptotically with the positive frequency 
WKB solutions fk as t ^ -oo. Thus |0) is the "out" 
vacuum in the parlance of [441 145j , and we assume that 
our system has been prepared in this vacuum. All expec- 
tation values unless otherwise indicated are with respect 
to this state. 
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B. The scaling limit 

In the classical context, we first solved the complete dy- 
namical problem for arbitrary initial condition and pro- 
tocol choice. We then established that all local physical 
quantities reduced to universal forms independent of the 
details of the protocol away from the critical point and 
the initial conditions in the KZ scaling limit. We repeat 
the same exercise here by using WKB methods to solve 
Eq. (|48| for a general TCP/CCP in m'^{t;T). 

Assume that the system is initially prepared in the vac- 
uum state. The solution to Eq. (48) is well-approximated 



by the positive oscillatory WKB solution fk at early 
times, and by a linear combination of the WKB solu- 
tions, fk and f^, at late times, as long as the frequency, 
1/Ofc, is much larger than the rate at which the fre- 
quency changes, \d\ogflk/dt\. This yields the condition, 
\dm~^/dt\ <C 1. Within the window |t| < r, we use the 
leading order expansion for m^(t; r) in t/r, 



m'^{t; t) = ni^ 



t 



1 



t 



(52) 



to conclude that the WKB solutions are valid as long as 
\t\ > t^ or \i\ > 1. Note that using Eq. Q is. 



t„ = T^+-/m^^ 



(53) 



Directly solving the mode equation using the expansion 
in Eq. (52) yields a solution valid in the region |i| ^ r 
or 1^1 ^ r/t^- Thus, the overlap in the ranges of the 
validity of the direct solution and the WKB one is: 



t. 



(54) 



We are now in a position to take the KZ scaling limit 
whereby r — oo and all quantities are measured in units 
of the KZ length and time. We first notice that the region 
of overlap in Eq. (54) diverges. The mode equation near 
t — also simplifies in scaled units. We see from Eq. ( 50 ) 
that fk carries dimensions of ^/t, so its scaling form is 



fk 



(55) 



where k = iKk as usual and Ik = tK- The scaled mode 
equation in the window |i| < r before any limits are taken 
is 

df2 



(56) 



where fi^ = 



m^{i) and 



l + ai- 



t„t 



02 



tKt 
T 



(57) 



In the KZ scaling limit, t^/r — and the corrections in 
the square bracket vanish, resulting in the simpler scaled 
mode equation. 



+ fc2 + (-£)" 



(58) 



The goodness of the approximation is parametrically con- 
trolled by the smallness of the parameter t^/r. Correc- 
tions to the diverging overlap in the ranges of the valid- 
ity of the solution to the scaled equation above and the 
WKB one is also controlled by the same small parame- 
ter tii/r. Thus, the solution to Eq. (58) picked out by 



starting with a positive frequency WKB solution at early 
times is the same as the one picked out by applying the 
early-time positive frequency condition directly to the 
solutions of the limiting equation, Eq. (58). This is the 



sense in which the mode functions in the scaling limit are 
universal. It is worth noting that these considerations are 
merely an elaboration of the standard arguments used to 
justify turning point formulas in standard WKB treat- 
ments of the time-independent Schrodinger equation in 
the geometric optics limit. 

Finally, let us note that the above considerations will 
not apply to the marginal ECP. For a given non-universal 
regularization of the protocol at small time (t < r) , the 
direct solution to the mode equation in the power-law 
regime {t ^ r) has no overlap with the positive frequency 
WKB solution at early times. As we discuss later, this 
implies that some information about the short time reg- 
ularization must enter the scaling limit. 



C. Quasi-particles, heat density and diagonal 
entropy 

As the Gaussian problem is integrable, quasi-particles 
are well-defined and infinitely long-lived. The quasi- 
particle number at momentum k, A4, is defined as the 
expectation value of the occupation of mode k and has 
the scaling form. 



^kii) 



(4Wak(i)) 
' ;\difk^ 



2\n;. 



(59) 



The excess energy density is given by q{t; r) = 
J d'^fc/(27r)''ilfeA4 and is easily seen to obey the scaling 
form conjectured in Eq. (25). A definition of the entropy 



density of the system at each instant of time is through 
the diagonal entropy density. The diagonal entropy den- 
sity, s, is the entropy density of the diagonal components 
of the density matrix p{t; r) in the many-body adiabatic 
basis. In the Gaussian problem, s is additive in the label 
k and this definition simplifies to. 



sit;T) 



d'^k 



i2n) 



d ^ 



Pmk.rnk l0g(pmk,mk) ■ 



Here Pmk.mk is the absolute value squared of the overlap 
between the time-evolved wave-function \ip(t)), and the 
nith excited state of the harmonic oscillator labeled by k. 
The integral only runs over half the volume in /c-space as 
the modes, k and — k, are coupled by the Hamiltonian. 
Pmk,mk Can be directly computed in the Schrodinger pic- 
ture because the time-evolved wavefunction is known in 
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terms of the mode functions. In the eigenbasis of the 
operator c^k defined in Eq. ( 46 1 , the wavefunction is a 
Gaussian liSl. 



bk\^) (X exp 




_^ .Aw 

\h(t)V "hit) 



up to normahzation and time-dependent phases. In the 
scahng hmit, fk{t) ^ V^iifki^) and s has the scahng 



form predicted in Eq. ( 25 1 . As all that is at issue is the 
scaling of the mode functions /fe(t), the scaling forms of 
correlation functions may also be easily verified. 

G^^{k,t-T) = |/fc(t)|2 ^ t^\J^(t)\^ ^ t^G^^ki) ■ (60) 



D. Linear protocol 

We now turn to a particular protocol, namely the linear 
quench, where the KZ time is = (t/toq)^^'^. For this 
case the mode equation can be solved in closed form to 
give 
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FIG. 5: G^if, for the linear TCP. The blue, red and green 
solid lines are respectively a,t t = 0.5, and —0.5. The green 
dashed line is the correlator if the system were in equilibrium 
at t = —0.5. Inset: Cnc vs t (solid) and ^/Ik (dashed). 



density from the high-wavenumber (fc 3> 
for t <C — I is finite only if d < 6: 



i) modes 



hit) = 



Bi(t- fc^) +i Ai(f- /c^) 



(61) 



where Ai and Bi are Airy functions of the first and second 
kind. The scaling function of the two-point function from 
Eq. jeOl) is: 



G(t,c/,ik,i) 



Ai" it 



Bi2(£-P) 



(62) 



As a check we may retrieve the equilibrium result when 
t — > — CXI holding or k/\/—t fixed. 



t/00(fc, t) 



1 

2% 



(63) 



The equilibrium result is not retrieved when i — >■ oo be- 
cause of the pathology of the inverted 0^ term when 
t > 0. Instead, grows exponentially with t. From 
G^^ we can calculate the non-equilibrium correlation 
length, 



i djiMiir+Biit)^) 

Ai(f)2 + Bi(£)2 



(64) 



whose scaling form £no is plotted in Fig. [sja). At large 
positive t, Cnc grows polynomially instead of exponen- 
tially with time, despite the instability of the adiabatic 
theory with an inverted 0^ term. As the adiabatic prob- 
lem is pathological for f > 0, we can only sensibly talk 
of quasiparticle occupations and thermodynamic quanti- 
ties as long as £ < 0. In Fig. [6] we plot A/'j.(t) and 5j,(t) 

for various values of k (the behavior of the excess energy 
density can be inferred from the number of quasiparticle 
excitations). The contribution to the total quasi-particle 



-Mi) 



1 



if d < 6. 



The total quasi-particle density, and consequently the en- 
ergy and entropy density, are ill-defined in the scaling 
limit within this integrable model when d > 6. Addition- 
ally, observe that diverges at £ = because the gap, 
17o, between the ground and excited state of the oscillator 
at fc = closes at t = 0. 



sdi) 




FIG. 6: Left: Quasiparticle number AAj, (£) for k = 0,0.25,0.5. 
Right: Diagonal entropy S^. for k = Q, 0.01, 0.1, 1 from top to 
bottom in the linear TCP. 



E. CCP - Quadratic protocol 

We now turn to a case where the adiabatic problem is 
well-defined for all t - the quadratic CCP. The KZ time 
here is = (t/toq)^/^. The scaling forms of the mode 
functions and the equal-time two point function are : 



hH) 



-1/4, 



i-V2e-Ti) 
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V2 



-|i?^i+^(-\/2e~Tt)|2 



(65) 



Dj^{z) is the parabolic cylinder function. 

As i — ?> ±00 holding /c^ or k/t fixed, we recover the 
equilibrium forms, 



1 

2^ 



(66) 



The retrieval of the equilibrium form as t — > 00 is by 
no means guaranteed. Recent work [25* suggests that 
this question is intimately tied to the de-phasing of the 
off-diagonal terms in the time-averaged density matrix 
in the instantaneous cigenbasis. If the off-diagonal terms 
do not de-phase, which would be the case for example 
for the quartic quench, then equilibrium behavior is not 
expected as i — )■ 00. 




0.0 0.5 



FIG. 7: Left: 6*^0 for the quadratic CCP. The blue, red 
and green solid lines are respectively at t — 0.5,0 and —0.5. 
The green dashed line is the correlator if the system were in 
equilibrium a.t t — —0.5. Right: Cno vs t (solid) and £,/Ik 
(dashed). 



This confirms the prediction in Eq. ( 18 1 for CCPs with 
no coarsening physics. For the reason discussed above, 
it also exhibits characteristic oscillations of period 1 /i at 
long times. 
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FIG. 8: Left: Quasiparticle number A/'fe(i) for k = 0,0.25,0.5. 
Right: Diagonal entropy Sf. for k — 0, 0.01, 0.1, 1 from top to 
bottom in the quadratic CCP. 



Finally, we plot A/'j,(t) and Sj^{t) for various values of 
k in Fig. [s] In each oscillator, the quasiparticles and 
entropy is essentially produced in the time- interval \t\ < 
1. At late times, both quantities settle to a constant 
dependent on k. For example, 



as t 



As was the case for the linear quench, the total quasipar- 
ticle and entropy density is finite only if d < 6. 



F. The Marginal ECP 



The Gaussian theory imposes further structure on cor- 
relators because each momentum mode evolves indepen- 
dently of the others. Consider the time-evolution of the 
wavefunction of the harmonic oscillator labelled by k in 
the Schrodinger picture. Expanding it in the eigenbasis 
of the Hamiltonian at time t and suppressing the label fc, 

W ^ AolOii)) + A,\l{i)) . . . . 

The time-evolution is adiabatic when t 3> 1 and each 
\Ai\ approaches a constant. The relative phase between 

Ai+i and Ai on the other hand grows as : J* Cl^{i')dt' ^ 
P. This is the origin of the oscillations of period 1/t at 
large positive t at each k, seen for example in G,j,^{k,t) 
inFig.IJa). 

The simple definition of ^ne in Sec. |III A4| has to be 
modified as Q^^(k,t) has multiple poles in the complex 
fc-plane for late times. The pole with the smallest imagi- 
nary part determines the decay constant (1/^ne) over the 
longest length scales in the two-point function in real- 
space. The non-equilibrium correlation length is plotted 
in Fig.[7][b). As £ -> ±00, the envelope of behaves as, 



For the relativistic Gaussian problem with v = 1/2 and 
z = 1, the mass in the marginal ECP vanishes as mQ^r/t) 
at late times t ^ t. toot is a dimensionless number 
that will play an important role in the following. Let us 
consider a general protocol with the twin requirements 
that it match this late time behavior and that it lead to 
adiabaticity at early times —00 < t ^ r. At late times, 
the system is led to the critical point with ^,^t ^ t. Our 
naive expectation is therefore that the mode function is of 
a universal scaling form, fk{t) = \/tf{kt), in the scaling 
limit t — )■ 00, fc — )■ with kt held fixed. 

This expectation is violated on two fronts even in the 
Gaussian problem. The first, as we previously discussed 



in Section V B is that the marginal ECP is by construc- 
tion unable to entire "forget" its early time regulariza- 
tion. More precisely, observe that our choice of pro- 
tocols implies that \dm~^/dt\ ^ 1 when t ^ t while 
\dm^'^/dt\ — l/(mor) when t ^ t. The latter re- 
sult shows that for small enough k, it is not possible 
to reach the power law regime while remaining adia- 
batic, in contrast to our discussion of the TCP and CCP 
cases; thus we should expect that some non-universal in- 
formation must make its way into the putative scaling 
regime. The second and more striking violation is that 
of scaling. When moT > 1/2 and \dm~^/dt\ is closer 
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to obeying the condition of adiabaticity, scaling is vio- 
lated "mildly" and all physical quantities are periodic 
functions of log(t). When niQT — 1/2, scahng is violated 
logarithmically. The most dramatic violation of scaling 
is when < mor < 1/2. Here the adiabaticity condi- 
tion is strongly broken and the scalar field acquires an 
anomalous dimension. 

Let us now describe how these features emerge in the 
long time, small momentum form of the mode functions. 
The mode equation for t ^ r: 



i2 



fk{t)=0. 



(67) 



Define A = 



and 



I iX, if m,oT > 1/2 
h{X) = {0, if mor = 1/2 
A, if mor < 1/2. 



(68) 



Setting fc = 0, Eq. (67) is solved as 



■^/^(*^% (A)'-'"+^(A)"™ 

foit) - { if mor 7^ 1/2 

y<(l + ulog(^)), if mor = 1/2 . 

The particular A-dependent pre-factors are chosen to sim- 
plify the solution at all k. io and u are non-universal 
constants in the equation above. They depend on the de- 
tailed pre-asymptotic form of the protocol and are fixed 
by the solution to the differential equation with the full 
time dependance of the mass. 

Let us now turn to fc ^ 0. For times greater than r, 
the mode equation is solved by the linear combination. 



fk{t) = Vt [cj{k) Jh(^){kt) + CY{k)Yj,(^){kt)\ . 



(69) 



Jv{x) and Yy{x) are the Bessel functions of order v of 
the first and the second kind. The functions cj{k) and 
cy(fc) are presently undetermined except for the require- 
ment of smoothness in k. However, if we examine the 
region i ^ r at small kt, we can obtain their leading 
order behavior by requiring that they match smoothly 
to the k = forms presented earlier. This localizes the 
non-universality to the constants u, tg discussed already. 
For mor < 1/2 and k smaller than some non- universal 
ko, this leads to the leading behavior: 

Cy(fc) = (kto)^ 

cj{k) = u{kto)-^ -^^^cotnX. (70) 

TT 

For mor > 1/2, we simply replace A by iX. For mor = 
1/2, we find 

cyik) = —u 




FIG. 9: g^,t,(kt) and J^{kt) vs kt when A = 0.1 and u/tg 
-2.5i. 



cj{k) = 1 - ulog 



kt^e^ 



(71) 



7 is the Euler-Mascheroni constant. We emphasize that 
the mode-function and consequently all physical quanti- 
ties are well-approximated by Eq. ( 69 ) with the forms of 



cj and cy above for t r and k ko- 

We now turn to the physical implications of these so- 
lutions for an arbitrary choice of (u, io) that satisfy the 
constraint imposed by the Wronskian condition Eq. ( |50[ ) . 
First consider the case < mor < 1/2, which is strongly 
non-adiabatic in the sense that \dm''^ /dt\ = l/(mor) is 
greater than 2 and may be large. The mode function in 
Eq. ( 69 ) can be re- written in the form 



g+ikt) 



-2A 



9-{kt) 



(72) 



and g- are linear combinations of the two Bessel func- 
tions and involve the constant u. Observe that in the 
limit of large t {t ^ to) with kt fixed, the ratio of the sec- 
ond to the first term in the expression above decreases as 
{t/to)~'^'^ [55]. Thus, in the scaling limit, fk{t) simplifies 
to 



h{t)^u't^l^+^f{kt)^u't^'^+^ 



Jxjkt) 
(fci)^ ■ 



(73) 



u' = u/tQ above. The scaling forms predicted in 
Sec. |IIIB] rfo not hold. Instead, the form above is the one 
expected when the field has an anomalous dimension 
A. The scaling form of G^^ with the modified dimension 
of (h is 



G^4.{k,t) ^ t'+'^g4.4.{kt) . 



(74) 



In real-space, this implies that the equal-time two-point 
correlator decays as l/x'^^^^^^ at fixed x/t. Analogously, 
the quasi-particle number Ni^{t) has the scaling form 
t^^N{kt). The scaling functions are: 



g^^{kt) = \u'\^p 
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Af{kt) = |u' 



, \{kt)dktf + iV{kty + {moTrf\' 



They are shown in Fig. [9j Three comments are in order. 
First, in the marginal ECP, ^ ~ i and the time-scale 
for a change in m is the same as ^ = ^t. We therefore 
expect that the non-equilibrium correlation length 
also grow linearly in time. This is indeed the case; the 
oscillations in Q^^{kt) ^ t in Fig.|9]are of order one period 
and indicate a peak at r ^ i in the real-space correlator 
G(t>(t>{r/t). Second, Nk{t) in the scaling limit grows as t^^ 
indicating the strong breaking of adiabaticity. Finally, 
all scaling functions are known up to a multiplicative 
constant (u') in the scaling limit. 




FIG. 10: G^^{k,t)/t and Nk{i) vs kt when \ — i plotted at 
three equally spaced time on the log(i) scale : t = 10 (red), 
(green) and 10 e^"'^'/^ (blue). The dashed line is 
the adiabatic response, (it, to) are chosen to be (\/3/2, 1). 



markedly from the adiabatic response, show as a dashed 
line in the plot. The period oscillations in the two-point 
correlator imply that ^ t- The oscillatory behavior 
of Nk [t] is an indicator that the system is closer to being 
in equilibrium than for the strongly non-adiabatic case 
< TOoT < 1/2. In a sense, the factors of t'l^l = e'l^l'°s* 
can be regarded as introducing only "logarithmic" mod- 
ifications of scaling. 

The case toot = 1/2 displays the properties of both the 
cases discussed above. The relations in Eq. ( 71 ) imply 
the form 

fk{t) = ^t[g+{kt) +\og{tm)g^{kt)] . 

At late times with kt fixed, the second term dominates 
the first and the mode function has the scaling form 

fk{t)^uVt\QgtJo{kt) 

All scaling forms are thus modified by pre-factors of pow- 
ers of \ogt. The logarithmic violation of scaling is sim- 
ilar to the weakly non-adiabatic case m^r > 1/2, while 
the modification of the dimension of (j) is similar to the 
strongly non-adiabatic case. The quasi-particle number 
here diverges logarithmically in time. 

Finally, let us comment on the ECP in the transverse 
field Ising model in (1-1-1) dimensions which is famously 
a model of free fermions. Here v = 1 and z = 1 so the 
relevant power is now a = 1. The explicit solution of 
the fermionic mode equations exhibits the analogs of our 
weakly non-adiabatic regime with corrections to scaling 
that are periodic in logi [5D]. Connecting this behavior 
to our Gaussian results discussed above by continuation 
in the number of dimensions is an interesting challenge 
for future work. 



Let us consider the case m^T > 1/2, which is weakly 
non-adiabatic in the sense that \dm~^ /dt\ = l/(moT) is 
less than 2 and may be small. A re-writing of the mode 
function informs us that the scaling is violated only by 
phases: 



g+{kt) 



g-{kt) 



(75) 



These phases affect other physical quantities in the scal- 
ing limit, but in relatively mild ways as compared to the 
previously seen factors of t^. For instance, at late times 
holding kt fixed, 

G^0(fc, t) ^ t g^^kt, e^l^l i°s(*/*o)) , (76) 

where we have expressed t*''*'! — e*l^l'°st order to em- 
phasize that this quantity is periodic in log t. The scaling 
function Afk is also periodic in \ogt. The two "almost- 
scahng" functions and N'kit) are plotted in Fig 



VF 



for three equally spaced values of log t within a decade. 
The late time behavior of Q^^ in the marginal ECP differs 



1. Analogies to dSd+i/CFTd 

We showed above that in the long time limit of the so- 
lution to the mode-equation in the strongly non-adiabatic 
case Eq. ( 72 1 , the scalar field acquires an anomalous di- 



mension A. If we instead took the opposite limit of small 

time (r <C t <C to), the second term is more important 
that the first and the anomalous dimension of the scalar 
field is —A. Thus, the effective scaling dimension of cj) in 
dimension d is 



A4 



1 



± A. 



(77) 



at short and long times respectively. Readers famil- 
iar with the {A)dSd+i/{C)FTd literature, in particular 
[171 HH] , will note a similarity between the result Eq. ( 74 ) 
and the d-dimensional field theory Green's function ob- 
tained in the presence of a double-trace deformation. We 
can make the analogy closer by noting that the modified 
scalar field 

d-l 
(fi = t 2 (j) 
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obeys an equation of motion which foUows from the La- 
grangian 

1 



C = ^^-dctg^,i{-g^-d^^d,^ - AfV') , (78) 



where g^i, is the metric of de Sitter space, dSd+i- 



dx' 



(79) 



The length scale L is arbitrary. The de Sitter mass must 
be given by 



d^ - 1 



(80) 



in order for the equation of motion from Eq. ( 79 ) to agree 



with Eq. (67). The dimensions Eq. (77) are closely re- 



lated to the usual ones in dS/CFT: 



A+ = :,± 



^2 - 1 



(81) 



The difference arises due to an important distinction be- 
tween the late-time Green's function Eq. (74 1 and the 
d-dimensional field theory two-point function. Examples 
of the latter could be computed in dSd+i in terms of early 
time (that is, small t) properties of mode functions as- 
sociated to the "out" vacuum: that is, mode functions 
which are purely positive frequency at late times. Our 
computation is the reverse of this, in the sense that we 
investigate late-time properties of mode functions asso- 
ciated with an "in" vacuum. Late time corresponds to 
the deep interior of dSd+i (more precisely, it is a corner 
of the global covering space far from the boundary at 
t = 0). Another difference (of lesser consequence) is that 
in the normal parlance of dS/CFT |47j, time flow would 
be reversed, so that what we call t = is the far future 
while t — > cxD is the far past. 

In the window fct <C 1, /fc(t) is fc- independent and equal 
to fo{t): 



large mpT^ implies large M'^L? . This in turn implies that 
the Compton wavelength of the massive scalar is much 
smaller than the Hubble scale of dSd+i, which is precisely 
the condition one needs in order to justify a geometric 
optics approximation. 



VI. CLOSING REMARKS 

Our primary purpose in this paper has been to sys- 
tematize the universal content of the KZ problem and to 
emphasize that all physical quantities give rise to uni- 
versal scaling functions that span the entire crossover 
from equilibrium at early times to the late time state. 
We have presented model computations that bear out 
this logic. Experiments directed towards observing this 
broader scaling picture would be highly desirable. 

An obvious challenge is to extend such computations 
to more physically realizable problems where the field 
theories are not as simple. We will present some results 
on non-trivial, but not physically realizable field theories, 
via the AdS/CFT correspondence elsewhere. Another 
obvious challenge is to formulate a renormalization group 
procedure that makes the universality manifest — beyond 
the case of stochastic classical models with the associated 
functional integral formalism discussed here. 

A byproduct of our analysis has been the identifica- 
tion of an especially interesting ECP which is able to 
produce anomalous dimensions already at the Gaussian 
level, through a mechanism similar to the way anomalous 
dimensions emerge in (A)dS/CFT. A deeper understand- 
ing of this phenomenon and its examination in interacting 
contexts is also a fit subject for further exploration. 



Appendix A: Scaling form of f{t\T) in the Gaussian 
theory with Model A dynamics 

To determine the generating functional in one higher 
dimension, we follow the three steps prescribed in |38l 



(82) 



where b^p, are known in terms of u, to, and in particu- 
lar have a definite ratio. This setup now bears a strong 
resemblance to double-trace operator deformations in 
AdS/CFT gHj. In AdS/CFT, the relativistic conformal 
symmetry of the boundary theory is broken by the multi- 
trace deformation. The presence of different powers of x 
in the real-space Green's function dependent on the en- 
ergy scale l/to signals a similar breakdown of scaling. 
The same can be said of the weakly non-adiabatic case, 
but comparisons with a boundary field theory are harder 
in this case because the field theory would have to be 
non-unitary, similar to violations of the Breitenlohner- 
Freedman bound |49] in AdS/CFT. However, the dSd+i 



formulation ( 78 ) does offer some further intuition regard- 
ing the weakly non-adiabatic case: Eq. (80 1 shows that 



• Define 0^ as the solution of the equation of motion 
m Eq. ml for a given noise history ( : J4>c, (k, t) — 
C(k, t). is a linear operator in the Gaussian the- 
ory. 

• Rewrite Z as j d(^P{C) exp{J(t)(^). The noise distri- 
bution, P(C), is Gaussian. 

• Recognize that the probability distribution for 0^ 
is related to the noise distribution as Prj,{(j)Q) = 
P{J(bc)detiJ). 

The generating functional of correlation functions of (f) so 
obtained is 



Z[J,t] = J V(j)det{J) e 



/ d^kdt {-2\J<p{k,t)\'^ + J(k,t)<p{-k,t)) 

(Al) 
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where 



J ^ g-^+k^ +ro{t;T). 



(A2) 



The functional integral is Gaussian, and the free energy 
density is expressed in terms of the structure factor in 
Eq. (pOl) as 



_fk_ 

(2^ 



log [G-,\k,t;T)] 



The challenge in identifying /na, even in the time- 
independent setting, lies in subtracting cut-off dependent 
terms from / that are analytic in S. Here, it involves 
subtracting the cut-off dependent equilibrium contribu- 
tion at the CP, /(O, cxd), and two terms that are linear 
and logarithmic in the cut-off. On taking the KZ scaling 
limit of the terms remaining, we confirm the scaling form 
in Eq. (pTl): 



/na(i;r)^^^(£) 



1 

6^ 



dk k^ 



k dQ^^ 
•J dk 



2 -t) 
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